Diamagnetic response of cylindrical normal metal — superconductor proximity 
structures with low concentration of scattering centers 
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We have investigated the diamagnetic response of composite NS proximity wires, consisting 
of a clean silver or copper coating, in good electrical contact to a superconducting niobium or 
tantalum core. The samples show strong induced diamagnetism in the normal layer, resulting in a 
nearly complete Meissner screening at low temperatures. The temperature dependence of the linear 
diamagnetic susceptibility data is successfully described by the quasiclassical Eilenberger theory 
including elastic scattering characterised by a mean free path I. Using the mean free path as the 
only fit parameter we found values of I in the range 0.1 — 1 of the normal metal layer thickness diss, 
which are in rough agreement with the ones obtained from residual resistivity measurements. The 
fits are satisfactory over the whole temperature range between 5 mK and 7 K for values of c£n varying 
between 1.6 /an and 30 fim. Although a finite mean free path is necessary to correctly describe the 
temperature dependence of the linear response diamagnetic susceptibility, the measured breakdown 
fields in the nonlinear regime follow the temperature and thickness dependence given by the clean 
limit theory. However, there is a discrepancy in the absolute values. We argue that in order to reach 
quantitative agreement one needs to take into account the mean free path from the fits of the linear 
response. [PACS numbers: 74.50. +r, 74.80.-g] 



I. INTRODUCTION 



A normal metal (N) in good electrical contact with 
a superconductor (S) exhibits superconducting proper- 
ties as the temperature is reduced. First experiments 
on the proximity effect were reported by R. Holm and 
W. MeissnerEl, who observed zero resistance between SNS 
pressed contacts. Since then, many .-investigations on 
proximity effects have beenxarried outn. Recently, it has 
received a revived interests. Particularly, experiments 
on the magnetic response have demonstrated nontrivial 
screening properties, showing .hvsteretic magnetic break- 
down at finite external fieldsETEl as well as a presently 
unexplained reentrant effect at low temperaturesQ. 

First theoretical studies on the proximity effect were 
carried out by Coopeiu. For a diffusive proximity system 
the diamagnetic susceptibility in the linear regime was 
investigated by the Orsay Group on Superconductivity 
in the framework of the Ginzburg-Landau theoryu. In 
this limit it was found that the induced diamagnetic sus- 
ceptibility in N depends on temperature approximately 
as xn oc 7 1-1 / 2 . |-This was confirmed experimentally by 
Oda and NaganoEil In spite of that, the magnetic proper- 
ties of proximity samples discussed in Refs. and o could 
not be explained by the Ginzburg-Landau theory, since 
they show a much stronger temperature dependence of 
the diamagnetic susceptibility. 

The clean limit, which is defined for the elastic mean 
free path I — ► oo, was first studied by ZaikinEH for a 
finite system of ideal geometry with the help of the qua- 
siclassical theory. He predicted that, since in this case 



the current -field relation is completely nonlocal, the cur- 
rent in N is spatially constant and depends on the vector 
potential integrated over the whole normal metal. As a 
consequence, the magnetic flux is screened linearly over 
the normal layer thickness and for T — + the sus- 
ceptibility can reach only 3/4 of that of a perfect dia- 
magnet. For T ^> Ta he found xn oc exp(— 2T/T\), 
with the Andreev temperature Ta = hv-p This 
mesoscopic temperature scale originates from the energy 
of the lowest Andreev bound stateEj. In this ballistic 
case all properties of the proximity effect are determined 
by these bound states, which are coherent superpositions 
of electron and hole waves between consecutive Andreev 
and normal reflections. 

Within the framework of the quasiclassical theory, 
Belzig et aiJifj studied recently the magnetic response of 
a proximity coupled NS sandwich in the two limits, clean 
and dirty. They obtained numerical solutions of the cor- 
responding equations for a wide range of temperature, 
magnetic field and layer thickness. Furthermore, they 
tried a fit of the susceptibility data of one AgNb specimen 
with a relatively small mean free path (specimen lAgNb 
in this paper) within their dirty limit results. The fit was 
only successful at very low temperatures, but could not 
reproduce the high temperature data. The diamagnetic 
screening of the dirty limit was too big as compared to 
the experiment. Moreover, the temperature dependence 
of the susceptibility of NS proximity specimens with the 
longest mean free paths could not be explained within 
the clean limit theory. This limit gives a completely dif- 
ferent temperature dependence than the one observed in 
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the experiments. In addition, data at low temperatures 
reach almost perfect diamagnetic screening, whereas in 
the clean limit it should not exceed 75 % of full screen- 
ing. To close the gap between these two limits the linear 
magnetic response for arbitrary impurity concentrations 



was theoretically studied in Ref. 14. We show in this 



paper that these results can satisfactorily describe the 
experimental data. 

The theory of the nonlinear reponse was addressed 
in Refs. ^,^3[ and lj| A proximity NS sandwich with 
a finite normal layer thickness g?n at sufficiently low 
temperatures undergoes a first order phase transition 
both in the clean and in the dirty limit. At a break- 
down field a jump in the magnetization occurs. In the 
clean limit the temperature dependence of the break- 
down field is H\y{T) cx exp(— T/TjCj. We show in this 
paper, that the clean limit result agrees well with our 
experiments on Hb (see also Ref. ^|), which follow the 
same temperature dependence in relatively clean speci- 
mens (l/cki w 0.4 — 0.8). Also the experimentally found 
absolute values of the breakdown fields show a 1/g?n _ 
dependence, whicb-is in qualitative agreement with the 
clean limit theorylia. However, this is not the case for 
the shape of the magnetization curves and, as mentioned 
before, for the temperature dependence of the linear sus- 
ceptibility. 

We have fabricated a new set of CuNb and AgNb spec- 
imens similar to the ones reported in Ref. ||. The spec- 
imens were produced with an optimized annealing pro- 
cedure in order to achieve very high mean free paths. 
The ratio of the thickness of the normal layer to the ra- 
dius of the superconductor of the composite wires, djsi/rs, 
was also varied to investigate comparable samples with 
normal layers reaching from almost flat to rather curved 
cylindrical geometry. Thus the influence of the NS ge- 
ometry on the proximity effect was investigated. In this 
paper, we discuss these new samples together with older 
onestj, covering a wide range of parameters, and fit their 
diamagnetic response with the help of the quasiclassi- 
cal theory in an intermediate impurity regime between 
the clean and dirty limit. Here, we will not consider the 
reentrant effect found by Mota and co-workersQ but only 
note that a recent, theoretical study has addressed this 
new phenomenon!!!! For the nonlinear response we com- 
pare the temperature dependent breakdown fields with 
the quasiclassical clean limit result. 

This paper is organized as follows. In Sec. || we de- 
scribe some theoretical aspects of the magnetic response. 
In Sec. Ill we describe the sample preparation and the 



measurement apparatus. In Sec. IV the experimental re- 
sults on the linear magnetic response are presented and 
the data are fitted using the theory from Sec. [H| The 
experimental results on the nonlinear response are ad- 
dressed in Sec. ^ Finally we draw some conclusions. 



II. THEORETICAL ASPECTS 

yilhin the quasiclassical theory of superconductivi- 
tyEZrtJ the magnetic response of a planar NS-proximity 
structure was investigated for arbitrary impurity concen- 
trations by Belzig et ai.E-3. Below we summarize some ba- 
sic results of this paper. The theoretical system consists 
of a semi-infinite superconductor in perfect contact with 
a normal metal of thickness d^ and specular reflecting 
outer surface. The pair potential A is taken to be con- 
stant in the superconductor and to be zero in the normal 
metal. The mean free path £ and the Fermi velocity vf 
are assumed to be the same throughout the system. 

In Ref. [l4| a variety of regimes are discussed, where the 
proximity effect is different from the previously studied 
clean and dirty limits. The regimes differ by the relative 
magnitudes of the thermal length = hvp/2Trk-BT, the 
mean free path £, and the thickness g?n, an additional rel- 
evant length scale, which has to be taken into account. 
In the case of a finite normal layer thickness g?n, the clean 
and the dirty limit are restricted to much smaller param- 
eter regions than previously believed. The dirty limit 
holds for £ <C £,T,d^s only, if also the mean free path 
£ <C £o- Here £o ~ vp/A is the coherence length in the 
superconductor. On the other hand, the regime £ < £t 
belongs to the ballistic regime, if £ > d^. 

To treat the screening problem, a general linear- 
response formula was derived which yields a non-local 
current-field relation in terms of the zero-field Green's 
functions. These functions are obtained from the solution 
of the Eilenbcrger equation including an impurity selfen- 
ergy. In general this has to be done numerically, since 
analytical expressions can be derived only in the clean 
and in the dirty limit. Finally, the Maxwell equations 
have to be solved. 

The chosen geometry is such that the NS interface lies 
in the y-z plane. Then the current in the normal metal 
has the form 



j y (x) = — I K(x,x')A y (x')dx' 



(1) 



where A y is the vector potential in transverse gauge, if 
the magnetic field is applied in z-direction. The kernel 
K(x,x') contains an exponential dependence on the dis- 
tance \x— x'\. The range of the kernel is given by £ in most 
of the temperature and impurity regimes. For £ > £,t and 
£ -C exp(2dN/£r) that is, at high temperatures, it is 
given by £t- It has thus a strong temperature depen- 
dence which has to be contrasted to the case of a bulk 
superconductor, where the range of the Pippard kernel is 
only weakly temperature dependent. 

The prefactor of the exponential in the kernel, which in 
general depends on coordinates, is related to the local su- 
pcrfluid density. It introduces an additional length scale 
in the problem, the field penetration depth X(x, T) = 
\wf(x,T), where An = (47re 2 ?i/m) _1 / 2 is a London like 
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length in the normal metal and f(x,T) is a function of 
temperature and position. The interplay between the 
range of the kernel and the superfluid density determines 
whether the nonlocal form of the current field relation is 
important. This has strong consequences on the screen- 
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FIG. 1. Numerical results of the screening fraction of an 
NS proximity slab for An = 0.003dN and different impurity 
regimes as functions of T/Ta- The clean limit is given by a 
solid line and an Usadel theory result for £ — 0.1 cZn is given 
by a dash-dotted line. For the intermediate impurity regime 
three cases are given: I = IO^n, I = g!n, and £ — 0.1g!n. 

In Fig. |l| we reproduce some numerical results0 of the 
screening fraction p/d-^ = — 4-7txn in a normal metal layer 
for An = 0.003c?n- Five curves as a function of tempera- 
ture, normalized to the Andreev temperature, are given 
for different mean free paths I. 

For the clean limit, £ — > oo, at zero temperature, the 
screening can only reach 3/4 of —1/Air. This is due to 
the infinite range of the kernel, or in other words, the 
complete non-locality of the current-field relation. At 
a temperature T ss 5Ta, the screening is exponentially 
reduced by thermal occupation of the Andreev levelso, 
and at higher temperatures screening is almost negligible. 

For I — lOd^ and £ — d^ screening is enhanced in 
comparison to the clean limit at low and high temper- 
atures. At low temperatures screening is more effective 
for c£n < & < oo, since the range of the kernel is given 
by finite I. At high temperatures, where the range of 
the kernel is £r, screening is more effective than in the 
clean limit but still nonlocal. For I = Q.ld^, already at 
T ~ Ta screening is considerably reduced, since the su- 
perfluid density is suppressed because the mean free path 
is smaller than d^. 

In the diffusive regime, it can be shown that for 
£ 2 < Anc?n, the nonlocal relation is reduced-,to the local 
current-field relation of the Usadel theoryc^l. In partic- 



ular, for the case An = 0.003<iN corresponding approx- 
imately to the experiments discussed in this paper, the 
condition £ <C c?n has to be met. Indeed, in Fig. |l| we il- 
lustrate the non-negligible difference in screening even 
for the case I — Q.ld-^, between the local Usadel re- 
sult and the nonlocal approach with elastic scattering 
described in Ref. [hi]. For T < 5T A the Usadel result 
exhibits less screening, whereas for T > 5Xa it shows a 
much more pronounced screening tail. 

After these recent results, covering a wide range be- 
tween the clean and the dirty limits, we have now the 
possibility to analyze the magnetic response of our rel- 
atively clean NS proximity specimens. To compare the 
experiments with the theory we have used the full form 
(0) valid for all impurity concentrations and solved the 
screening problem for each specimen numerically. 

In these theoretical considerations the effect of a rough 
normal metal-vacuum boundary has been neglected com- 
pletely. On the other hand, at high enough temperatures 
T ^> Ta, where the induced diamagnetism strongly de- 
creases, the influence of surface scattering does not play 
an important role, since in this temperature regime the 
superfluid density at the outer boundary is exponentially 
suppressed. The screening in this regime can be only due 
to bulk scattering centers and is insensitive to the quality 
of the surface. 

The quasiclassical theory for the nonlinear response in 
the-clean limit was discussed in detail by Fauchere et 
allB. They derived analytical expressions for the tem- 
perature dependent breakdown field H\,(T) of a clean 
normal-metal slab of finite thickness in proximity with 
a bulk superconductor in the two limits T — > and 
T ^S> Ta- The breakdown field H\, was determined with 
a Maxwell construction from the magnetization at the 
bistable regime, where two different values of the free 
energy coexist, characterizing a diamagnetic and a field 
penetration phase. The spinodals of that bistable regime 
have already been determined in the clean and in the 
dirty limit by Belzig et alE^I numerically. They repre- 
sent the boundaries of superheating and supercooling in 
the first order phase transition. The thermodynamical 
magnetic breakdown field lies in between the spinodals. 
It was determined by the intersection of the two asymp- 
totics of the free energy. ._. 

For T — > the breakdown field saturates att£l 



«■(»'" stt- 

An UN 

whereas for T Ta it decays asEl 



H h (T) 



V2 $, 



7T An^N 



(2) 



(3) 



The temperature dependence is a simple exponential 
with the exponent c?n/^t = T/Ta- The amplitude of 
the breakdown field was found to scale inversely propor- 
tional to the normal layer thickness. Both features are in 
agreement with previous experimental results on H\> in 
relatively clean AgNb specimens^. 
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In the H-T phase diagram a critical temperature T cr it 
was found, below which the first order phase transition 
is observable. Above T cr ; t a continuous and reversible 
crossover between the diamagnetic and field penetration 
regime occurs. 

Since the breakdown field depends only on thermody- 
namical considerations, it is expected to be less sensitive 
to the nonlocality of the current-field relation than the 
linear susceptibility. A cylindrical geometry and a barrier 
at the NS interface is expected to change the breakdown 
field only quantitativelylij. On the other hand, the in- 
fluence of scattering centere is reflected in the shape of 
the magnetization curveal3. At present, no theoretical 
results for the breakdown field or magnetization curves 
exist for intermediate impurity concentrations. 



III. SAMPLE PREPARATION AND 
EXPERIMENTAL SETUP 

The samples we investigated are bundles of cylindrical 
wires with a superconducting core of niobium or tanta- 
lum concentrically embedded in a normal metal matrix 
of copper or silver. The normal metal starting materi- 
als are of purity 4N to 6N with negligible contribution 
of magnetic impurities. The purity of the niobium and 
tantalum starting materials was the highest accessible 
(RRR « 300), resulting in a good ductibility. We as- 
sembled the samples by placing in a close fit, a rod of 
the S material inside a hollow cylinder of the N material, 
with typical outer diameter of about 20 mm. The cleanli- 
ness of the internal N tube surface and the S rod surface 
is very important. The surfaces of the initial metal pieces 
were cleaned in acid and mechanically smoothed with a 
special tool in order to provide a good metallic contact. 
All this was done in an argon atmosphere. The NS cylin- 
der was then embedded in a Cu mantle to protect the 
sample itself from contamination. After assembly, the 
diameter was reduced mechanically by several steps of 
swagging and co-drawingEj down to a total diameter of 
several hundred microns. 

After etching away the Cu protection, the diameter 
of the NS specimens was further reduced by co-drawing 
the wires through several diamond dies to final values 
between 15 /im and 190 /im. The extreme size reduction 
(by factors of up to 1000) by co-drawing resulted in a 
great enhancement of the quality of the NS interface for 
electronic transmission. The use of high purity starting 
materials guaranteed a high electronic mean free path. 
Some of the samples were annealed after the last draw- 
ing in an atmosphere of argon at a temperature around 
650° C for about one hour to remove the effect of cold 
working in the normal metal. For that purpose the wire 
of selected size was wound on a small Ag plate. As a con- 
sequence of the optimized annealing procedure, values of 
^ f^N could be achieved. The mean free paths were 
obtained from resistivity measurements along the wires, 



performed at a temperature T = 10 K with a four-point 
method. The thickness of the normal metal layer as well 
as the surface quality of the wires were determined with 
the help of SEM micrographs. Several samples were pre- 
pared in order to cover a large range of parameter space. 
We fabricated AgNb samples with a ratio d-^/rg of 0.4 
or 1. The layer thicknesses ranged from c?n = 2.8 /im 
to c?n = 28 fim with measured mean free paths ranging 
from ^n/^n = 0.12 to £n/^n = 1-3. Moreover, we pre- 
pared CuNb and CuTa samples, with varying ratios of 
N layer thickness to S core radius, and mean free path 
£. This gave us the opportunity to apply the theory de- 
scribed in Sec. || in a wide range of samples. The values 
of the sample parameters are listed in Table Q. 

As a last step, the wires were glued with GE 7031 var- 
nish and thus electrically insulated. After drying, they 
were removed from the silver support, cut to a length 
of typically 3 mm to 5 mm and rolled together forming a 
bundle of 200 to 800 wires. The wire bundle was then 
placed directly inside the mixing chamber of a dilution 
refrigerator in contact with the liquid 3 He- 4 He solution, 
mounted parallel to the coil axis of the magnetic ac and 
dc field. 

The ac magnetic susceptibility was measured at tem- 
peratures between 5 mK and 7 K using an rf-SQUID sen- 
sor, inductively coupled via a dc-flux transformer to the 
sample. A certain fraction of the voltage applied to the 

TABLE I. Table of sample parameters. The measured 
value ^n/o!n and the fitting parameter t/dn are given. In 
the last column the saturated superheated field at the lowest 
temperature is shown. 



sample 


d N 


ratio 


Tan ii 


£n 


T A 


An 


£ N /d N 


£/d N 


tj sat 




[/im] 


d N /rs 


[°C] 


[/im] 


[mK] 


H 






[Oe] 


lAgNb 


14.5 


0.4 


not 


1.8 


120 


22 


0.12 


0.15 


1.62 


2CuNb 


15.3 


1.0 


not 


0.9 


120 


18.3 


0.06 


0.025 


0.42 


3AgNb 


5.5 


0.4 


700 


4.3 


310 


22 


0.78 


0.65 


15.2 


4AgNb 


6.8 


0.4 


700 


5.1 


250 


22 


0.75 


0.75 


12.9 


5AgNb 


3.3 


0.4 


800 


1.6 


520 


22 


0.49 


0.55 


24.7 


6AgNb 


28.0 


0.4 


550 


9.0 


61 


22 


0.32 


0.25 


2.07 


7CuNb 


5.3 


1.0 


not 


0.9 


360 


18.3 


0.18 


0.08 


5.9 


12 AgNb 


3.6 


0.4 


not 


1.8 


480 


22 


0.5 


0.26 


15.4 


14CuTa 


5.0 


1.0 


not 


1.0 


380 


18.3 


0.2 


0.07 


7.0 


16AgNb 


3.3 


0.4 


800 


1.7 


520 


22 


0.5 


0.35 


24.5 


19CuTa 


3.8 


1.0 


600 


2.5 


510 


18.3 


0.67 


0.13 


12.9 


20AgNb 


2.7 


0.4 


not 


1.4 


630 


22 


0.53 


0.24 


20.2 


21 AgNb 


3.6 


0.4 


800 


1.0 


480 


22 


0.28 


0.75 


15.8 


23CuNb 


5.7 


0.2 


not 


1.4 


340 


18.3 


0.24 


0.11 


8.4 


24CuNb 


2.5 


0.2 


not 


0.9 


760 


18.3 


0.35 


0.15 


23.9 


25CuNb 


1.6 


0.2 


not 


0.5 


1200 


18.3 


0.27 


0.19 


40.5 


28AgNb 


14.5 


0.4 


700 


5.6 


120 


22 


0.39 


0.45 


4.89 


33AgNb 


7.7 


1.0 


700 


5.3 


220 


22 


0.7 


0.35 


10.2 


34AgNb 


3.1 


0.4 


700 


4.0 


540 


22 


1.27 


0.55 


26.2 


35AgNb 


2.8 


0.4 


670 


3.6 


620 


22 


1.3 


0.55 


30.8 


36AgNb 


5.1 


1.0 


670 


5.4 


330 


22 


1.05 


0.4 


14.1 


37AgNb 


4.2 


0.4 


670 


3.5 


410 


22 


0.84 


0.45 


18.9 


38AgNb 


3.8 


0.4 


not 


1.2 


450 


22 


0.32 


0.18 


10.6 
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primary ac-coil was also mutually fed into the flux trans- 
former loop, using the SQUID as a null detector. This al- 
lowed us to measure the susceptibility of the samples with 
a relative precision of about 10 . Typical ac amplitudes 
H ac are between 0.06 mOe and 33 mOe, and frequencies 
16 Hz, 32 Hz, 80 Hz, and 160 Hz. The temperature was 
measured through the Curie-type magnetic susceptibil- 
ity of the paramagnetic salt CMN (cerium magnesium 
nitrate), which was calibrated with two Ge resistors. The 
methods and whole apparatus are explained in more de- 
tail in Ref. %. 



IV. LINEAR MAGNETIC RESPONSE 

A. Experimental results 

In the following we report on the temperature depen- 
dent ac-susceptibility without external dc-field, which 
we will address as linear response. We show here the 
results for some typical samples, which represent differ- 
ent behaviors and sample parameters, namely annealed 
silver-niobium samples with a large measured mean free 
path £n as well as not annealed ones with smaller 
Some of them have a ratio d-^/rg of 1, the normal layer 
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FIG. 3. Susceptibility data (o) and numerical fits (solid 
lines) of silver niobium samples as functions of T/Ta- Fitting 
parameters (a) i/dn = 0.25 and (b) £/d^ = 0.55. 
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FIG. 2. Susceptibility data (o) and numerical fits (solid 
lines) of silver niobium samples as functions of T/Ta- (a) 
Fitting parameter i/dn = 0.55. The dashed line denotes the 
clean limit and the dotted line denotes the i/dti = 0.4 curve 
for comparison, (b) Fitting parameter l/dvi = 0.25. 
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FIG. 4. Susceptibility data (o) and numerical fits (solid 
lines) of silver niobium samples as functions of T/Ta- Fitting 
parameters (a) l/da = 0.75 and (b) l/d^ — 0.55. The insets 
are zooms in the low temperature region for both samples. 
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FIG. 5. Susceptibility data (o) and numerical fits (solid 
lines) of silver niobium samples as functions of T/Ta- (a) 
ratio dfi/rs = 0.4 and fitting parameter l/d-^ = 0.75. (b) 
ratio d^/rs = 1 and fitting parameter £/dw = 0.4. The fit for 
sample 33AgNb shows a bigger deviation from the data for 
5 < T/T A < 15 than the fit for sample 4AgNb. 

having a rather strong curvature, while others have a ra- 
tio 0.4 and therefore an almost flat normal layer. Also 
shown are not annealed copper-niobium samples with ra- 
tios of 1 and 0.2 and a relatively smaller mean free path 
£n as well as two copper-tantalum samples with a ratio 
d-N/rs = 1. 

The data of the diamagnetic screening fraction 
p{T)Jd-N of the N layer are shown in Figs, g |, § |, 
and | as a function of temperature in units of the An- 
dreev temperature Ta for each sample. The rather strong 
development of the induced superconductivity in the nor- 
mal metal right below the transition temperature T c of 
the superconductor can be observed quite clearly. 

From the susceptibility x(T) measured in arbitrary 
units the screening fraction was obtained as 

p(T)/d N = r s /d N ■ [(1 + xn(T)/A Xs (T c ))^ 2 - 1] . 

Here 

XN (T)/A X s(T c ) = ( X (T) ~ XoVUo - Xoo) 

is the fraction of the temperature dependent susceptibil- 
ity of the normal metal with respect to the total diamag- 
netic transition of the superconductor at its critical tem- 
perature T c . The exact height of both A%s and Xn(T) 
depends on the value of the susceptibility xo right below 
the transition of the superconductor. Since that transi- 
tion for the niobium samples is outside of our accessible 




FIG. 6. Susceptibility data (o) and numerical fits (solid 
lines) of (a) the copper niobium sample 24CuNb, fitting pa- 
rameter £/dn = 0.15 and (b) the copper tantalum sam- 
ple 19CuTa, fitting parameter l/d^ — 0.13 as functions of 
TIT a- 

temperature range, values at the highest temperatures 
were used for extrapolation. The susceptibility above the 
superconducting transition Xoa was obtained from back- 
ground measurements. Through the ratio of the radius of 
the superconductor rg to the normal layer thickness g?n 
the normal layer curvature of the composite cylindrical 
sample is taken into account. 

In a local picture as given by the Orsay Group on 
SuperconductivityO the screening fraction p/d^ repre- 
sents the screening length normalized to the normal layer 
thickness, with p the thickness of the part of the normal 
metal, out of which the magnetic flux is screened. For our 
samples that picture is not valid, since the current-field 
relation is nonlocal. Here the screening fraction p/d^ 
represents the susceptibility of N, as related to a model 
of a one dimensional system. 



B. Discussion of the fits 



The temperature dependence of the diamagnetic sus- 
ceptibility of these samples is neither accounted for by 
the clean limit nor the dirty limit theory. The values of 
the measured mean free path £n indicate that the sam- 
ples are in an intermediate impurity regime. They are 
not diffusive enough for the dirty limit, because £n is of 
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the order of c?n for most of the samples. However, since 
they have a small density of scattering centers, they are 
not in the clean limit, either. From here on we will ad- 
dress them as relatively clean. 

We fitted the experimental data of the screening frac- 
tion with the theoretical results obtained for a one- 
dimensional geometry, as described in Sec. ||. The fits 
are shown in Figs, g, ||, ||, ||, and [|. 

The parameters c?n and An entering the numerical 
calculations were obtained from measurements and well 
known material constants. The only fitting parameter 
was then the mean free path £, which enters into the the- 
ory through the fraction £/d^- The fits were performed 
giving the experimental data at an intermediate temper- 
ature regime the most weight, where the increase of the 
susceptibility was the steepest. To illustrate the sensi- 
tivity of the fits to the value of the fit parameter £/d-^, 
in Fig. ||(a), we give in addition to the fitted curve with 
i/d-n = 0.55, the curve for £/d N = 0.4. 

The curves obtained from the newly developed ap- 
proach with arbitrary scattering center concentration fit 
the experimental data well over the whole temperature 
range. The corresponding clean limit curve does not fit 
the data, as shown for comparison in Fig. ^|(a). 

The values of the fitting parameter l/d^ reproduce the 
measured values £n/^n rather well. It has to be em- 
phasized, that the fits are rather good, if one takes into 
account the difference in geometry between experiment 
and theory, the neglect of the boundary roughness in the 
theory and other imperfections inevitably present in the 
samples. This means, that the quasiclassical theory of 
the proximity effect with a finite mean free path par.aHie- 
ter I due to a low concentration of elastic scattereraia, is 
now able to explain the linear susceptibility data of our 
relatively clean NS specimens. A list of all the samples 
shown here is given in Table [j]. Additionally more sam- 
ples are listed together with the mean free path I which 
gave the best fit. 

For the silver-niobium samples with a measured mean 
free path £n/^n ~ 1 the fitting parameter £/d^ lies be- 
tween 0.4 and 0.8. The measured mean free path of AgNb 
samples with ^n/<^n < 0.3 is also reproduced quite well 
in the fits. A typical example is the silver-niobium sam- 
ple 6AgNb shown in Fig. |^(b). It is suprising, that for 
most of the samples the agreement between the only fit- 
ting parameter £/d^ and the measured £n/^n is so good. 
However, the parameter £ has to be viewed as an effective 
mean free path, which contains the scattering from bulk 
impurities as well as from the surface. The measured 
mean free path £n is also an effective quantity, which 
contains surface effects, but it is yielded from transport 
measurements along the axis of the wires. There the 
surface is expected to influence the mean free path in a 
different way, because the geometry of the relevant tra- 
jectories differs from the ones in the proximity case. 

In the following, we describe as an example, the quality 
of the fits for the silver-niobium samples. These samples 
show the best agreement between the measured magnetic 



susceptibility and the theory. The results for some of 
them are illustrated in Figs. ||, |[ ||, and ||. In general, 
the fits given by the solid lines, describe the experimen- 
tal data well. In particular, for some of these samples 
we observe considerable deviations at the lowest temper- 
atures (T < Ta) where the theoretical curves saturate as 
shown in Fig. [I]. Moreover, the deviations are different for 
nominally similar samples. For example, sample 2lAgNb 
(inset of Fig. ||(a)) show higher x~n values than the the- 
ory, while sample 5 AgNb (inset of Fig. f|(b)) shows also 
higher values and in addition a strong reentrance of \n 
at the lowest temperatures. These apparently contradic- 
tory results are possibly related to the surface quality. At 
this point, it is not clear how a non-ideal reflecting sur- 
face affects the susceptibility, neither experimentally nor 
theoretically. In addition, silver-niobium samples with 
roughly the same c?n and £ show strongly different levels 
of reentrance. The origin of these different behaviors will 
be subject of further investigations. 

The fit shown in Fig. ||(b) is rather poor, which could 
be partially explained by the geometrical non-uniformity 
of sample 35AgNb, showing effectively two different val- 
ues of c?n (see also Fig. ^ and Sec. |v|). 

For T > Ta the fits of the AgNb specimens are quite 
good. For all the specimens, in a medium temperature 
regime, the theoretical susceptibility is a bit too high 
with respect to the measured one. This behavior can be 
observed more clearly for the samples in Fig. [g. Their 
normal layer thickness with values of d^ — 6.8 fim and 
c?n = 7.7 /im as well as the measured mean free path with 
values of = 5.1 fim and = 5.3 /jm are approximately 
the same. The only difference between the two samples 
is the ratio d-^/rs, which for sample 4AgNb is 0.4 and for 
sample 33 AgNb is 1. As expected, the sample 33 AgNb 
with the bigger ratio and therefore with the stronger cur- 
vature of the normal layer shows the more pronounced 
deviation in the medium temperature regime. 

The susceptibility of the copper-niobium samples is 
also described by the theory rather well. A typical ex- 
ample is shown in Fig. |^(a). For similar samples the 
fitted values of £/d-N are very small, lying between 0.02 
and 0.2 (see Table Q). In contrast to the silver-niobium 
samples with very small mean free paths, for the copper- 
niobium samples the value of £/d^ is only about half of 
the measured value of €n/^n- The CuNb samples show 
similar deviations in the different temperature regimes as 
the silver-niobium samples. 

The poorest agreement between experiment and the- 
ory is met for the copper-tantalum samples. An example 
is shown in Fig. ||(b). As the temperature is increased 
the theoretical curve starts to show much more screen- 
ing than the experiment. The deviation starts to show 
up, where the induced diamagnetism is about half re- 
duced. At higher temperatures the difference between 
the two curves is even larger. Moreover, the gained fit- 
ting parameter £/d^ is much smaller than the measured 
£n/<^n- Interestingly, the two copper-tantalum samples 
investigated (14CuTa and 19CuTa) show the same strong 
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disagreement between £ and £n- This suggests that some 
other effects should be considered for the combination of 
these two materials. 



V. NONLINEAR RESPONSE MEASUREMENTS 

We have investigated for our NS proximity samples 
the nonlinear magnetic response. With the experimen- 
tal setup described in Ref. ^| ac-susceptibility measure- 
ments at constant temperature as a function of a dc 
magnetic field were performed, as well as isothermal de- 
magnetization curves. 



(a) 
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FIG. 7. (a) Nonlinear ac susceptibility as a function of mag- 
netic field H. (b) Isothermal dc magnetization curve. 

Fig. ^ shows the nonlinear susceptibility and the 
dc-magnctization curve for one of the cleanest silver- 
niobium samples at T — 8 mK. The magnetic break- 
down of the induced superconductivity occurs at a mag- 
netic field iib ~ 25 Oe. Below the breakdown field 
the magnetic flux is totally screened, and above, it en- 
ters into the normal metal layer. The first order phase 
transition shows the features of superheating at increas- 
ing field and of supercooling at decreasing field. That 
kind of hysteretic behavior was already observed and 
discussed in the .dirty limit by the Orsay Group on 
Superconductivity^. This particular sample shows a 
small second transition, which is due to the niobium core 
not being perfectly centered resulting in effectively two 
slightly different normal layer thicknesses. We notice that 
the magnetic breakdown transition is not infinitely sharp. 



In our samples consisting of several hundred wires the su- 
perheated and supercooled transitions are triggered for 
each individual wire at slightly different fields, leading to 
a statistical broadening of the breakdown jumps. 

From the isothermal XN(-£0 - curves we have deter- 
mined the supercooled and superheated fields as the fields 
at the middle of each transition. In this way values of H sc 
and i? s h at different temperatures were obtained. 

For our relatively clean silver-niobium samples at tem- 
peratures higher than the Andreev temperature Ta °c 
wf/^N) the experimental breakdown fields follow the ex- 
ponential dependence -ffb(T) oc 1/c£n exp(— g?n/£n(T)). 
The experimentally found values of £n (T) reproduce the 
theoretical clean limit coherence length in silver £t = 
hv-F fluk^T = 1.69/im/T(K) within a few percent. For 
the sample 35AgNb this is illustrated in Fig. ||, where 
the experimental supercooled and superheated break- 
down fields as a function of temperature are given. The 
theoretical clean limit breakdown field of a one dimen- 
sional NS slab for T 3> Ta is shown as a solid line. A 
factor of about 0.3 has been used to shift down the the- 
oretical curve to fit the experiment. We notice that the 
temperature dependence of the breakdown field agrees 
qualitatively with the clean limit resultiLa, obtained for a 
one dimensional NS slab assuming ideal boundary condi- 
tions, e.g. ideal normal electron transmission at the NS 
interface. Quantitatively, the measured breakdown fields 
are a factor of about 0.3 smaller than the high tempera- 
ture theoretical curve in Fig. ^. 

This discrepancy could have different reasons. First, in 
our samples normal reflections at the NS interface may 
happen due to the mismatch of the Fermi velocities be- 
tween N and S and impurities or interdiffusion. Second, 
our samples have a cylindrical geometry, which weakens 
the proximity effect with respect to the one dimensional 
flat geometry considered in the theory. And third, the 
influence of impurity scattering on the degree of nonlo- 
cality is expected to lead to quantitative deviations of 
the breakdown field data from the clean limit theory. 

In Fig. [)] the saturated superheated fields of the silver- 
niobium samples are plotted versus the inverse normal 
layer thickness 1/cZn- The samples with £n/^n ~ 1 and 
d^/rs = 0.4 show a clear l/^N-dependence of the break- 
down field, in agreement with the clean limit theorycJ. 
On the other hand, the absolute values of the saturated 
breakdown fields are about a factor of 1.7 lower than in 
the clean limit theory. All the other samples with a ra- 
tio dw/rg = 1 or with a lower mean free path show 
smaller breakdown fields. As already mentioned, for the 
copper samples the influence of a reduced normal elec- 
tron transmission coefficient at the NS boundary could 
be especially strong leading to a further suppression of 
the breakdown field. 

The silver-niobium samples with £n/^n < 0.3 and 
the copper samples also show a temperature dependence 
Hb(T) oc exp(— g?n/£n(T)), but for these samples exper- 
imentally we found-i^N(T) = p- £r, where the prefactor p 
is about 0.3 to 0.93. The breakdown field data of these 
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samples can no longer be explained by the clean limit 
theory, since the impurity scattering is too strong. 
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FIG. 8. Measured supercooled (o) and superheated (•) 
breakdown fields of the annealed silver-niobium sam- 
ple 35AgNb with d]si — 2.8 /im. The line represents the clean 
limit theory of a one dimensional NS slab with a correction 
factor 0.3. 
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9. 1/c£n -dependence of the saturated superheated 



field for silver-niobium samples with ^n/o!n ~ 1 (•)• Slope 
of solid line: 88 Oe^im. The dashed line represents the clean 
limit result, with slope: 150 Oe/im. Also shown are less clean 
silver-niobium samples (In/^n < 0.3,o) and copper samples 
(x), which we do not consider in the fit. 

Although the temperature dependence of the break- 
down field is the same as given by the clean limit, the 
full form of the magnetization curves M{H) of the clean- 
est AgNb samples can not be explained within this limit. 
The slope below the transition does not correspond to 
the —3/4 of a perfect diamagnet and additionally, the 
observed finite screening tail above the transition (see 
Fig. 0(b)) does not appear in the clean limit. Clearly, 
agreement between the quasiclassical theory and the 
magnetization data could be achieved by taking into ac- 
count the finite mean free path t. (Note the similarity 
between the magnetization curve in Fig. 0(b) and the the- 
oretical dirty limit magnetization curves from Ref. [l^, 
which have the right slope as well as a finite screening 



tail, but give the wrong field scale.) Theoretical work on 
M(H) for arbitrary low impurity concentrations is still 
needed. 



VI. CONCLUSIONS 

We investigated the magnetic response of NS proxim- 
ity layers of high purity, and compared the experimental 
data with the quasiclassical theory. A large number of 
samples containing combinations of Ag and Cu with Nb 
and Ta has been fabricated and measured. All the sam- 
ples show an induced diamagnetic susceptibility close to 
— 1/47T for temperatures T < Ta. The mean free path £n 
obtained from resistivity measurements varied between 
O.ldN and 1.3dN. With the results of Belzig et al. which 
are based on the quasiclassical theory including arbitrary 
impurity concentrations, we were able to reproduce the 
experimental data quite well, using a mean free path I 
as the only fitting parameter. The mean free path I 
obtained in this way agrees within a factor of 2 with 
the mean free path 1^ determined by resistivity measure- 
ments. This good agreement shows that the linear dia- 
magnetic response of a proximity system with arbitrary 
impurity concentration is very well described with the 
present nonlocal approach. However, the very low tem- 
perature behavior, where some deviations occur and in 
addition an unexplained reentrance of the susceptibility 
appears, is still not understood. 

Magnetization measurements show a first order transi- 
tion at a breakdown field from a state with almost com- 
plete flux expulsion to a state with weak screening. The 
temperature dependence of the breakdown field for the 
cleanest specimens is in accordance with the clean limit. 
Nevertheless, the complete magnetization curve, in agree- 
ment with the linear susceptibility data, reflects the finite 
mean free path of the samples. 

With this work we show that the proximity theory, 
based on the quasiclassical approximation including the 
full nonlocal current response, can successfully describe 
the magnetic response of very clean samples. The ex- 
perimentally unavoidable low level of impurities can be 
accounted for with a mean free path £ which is of the 
order of the sample size. 
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